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Abstract. 

We study a model of dilute oriented loops on the square lattice, where each loop is 
compatible with a fixed, alternating orientation of the lattice edges. This implies that 
loop strands are not allowed to go straight at vertices, and results in an enhancement of 
the usual 0(n) symmetry to U(n). The corresponding transfer matrix acts on a number of 
representations (standard modules) that grows exponentially with the system size. We derive 
their dimension and those of the centraliser by both combinatorial and algebraic techniques. 
A mapping onto a field theory permits us to identify the conformal field theory governing 
the critical range, n < 1. We establish the phase diagram and the critical exponents of 
low-energy excitations. For generic n, there is a critical line in the universality class of the 
dilute 0(2n) model, terminating in an SU(n -I- 1) point. The case n = 1 maps onto the 
critical line of the six-vertex model, along which exponents vary continuously. 


1. Introduction 

Loop models and their critical universality classes play a major role in several areas of 
theoretical physics, such as the study of geometrical statistical models, quantum integrable 
models, conformal held theory [^, or quantum information theory |^. The simplest—and 
most studied—of all the loop models is conveniently described by starting from an oriented 
square lattice, represented in hgure [T] 
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Figure 1: Left panel: oriented square lattice used for the definition of the various loop models 
in this paper. Right panel: conhguration of completely packed, self-avoiding and mutually 
avoiding loops on this oriented lattice. 


Every edge of the lattice is then supposed to be occupied by a loop segment (a monomer), 
while vertices can be split in two possible ways compatible with the link orientations. This 
gives rise to conhgurations of completely packed, self-avoiding and mutually avoiding loops 
(see hgure[^. Finally, a loop fugacity n is introduced, giving to each conhguration a weight 
ri'^ioop®, where iVioops is the number of loops. 

It is well-known that this model has an underlying U(?7,) symmetry when n is an 
integer [^|^. This can be seen in two ways. The hrst stems from calculating the partition 
function of this model by a transfer matrix approach, with the x-axis representing the ‘space’ 
and the y-axis the ‘imaginary time’ direction. The transfer matrix is then a product of 
elementary operators at every vertex which produce the two possible splits that respect the 
lattice orientation. Graphically, each vertex can be represented as follows: 

X ,X ■ X -X 

or, symbolically, 

t = I + xe , (1) 

where the e, are generators of the Temperley Lieb (TL) algebra, acting on strands number i 
and z -|- 1. They satisfy a well-known set of multiplication rules that can be read from their 
geometrical definition [^|^. In particular. 


(g)^ = 


0 

r\ 


= n 


KJ 

r\ 


= noi, 


( 2 ) 


123 i N 123 i N 

so as to give every loop a weight n. The parity of the number of lattice sites N will play 
some role in the following. In the case of periodic boundary conditions horizontally, the 
consistency of the orientations of edges (see figure requires N to be even. With free 
transverse boundary conditions N can have any parity, but in most of the paper (unless 
explicitly stating the contrary) we shall nevertheless set N = 2L also in that case. 
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The TL algebra admits a simple Hilbert space realisation where links oriented to the 
left (resp. the right) on a horizontal line carry a fundamental n (resp. conjugate fundamental 
h) representation of U(n)||] The generators Cj can then be written as Cj = nP, where P is 
the projector onto the identity in the product n ® h (resp. fi ® n). The transfer matrix (or, 
in the anisotropic limit, the Hamiltonian) then acquires a U(n) symmetry. 

General arguments mapping the spin chain Hamiltonian to a sigma model (see, e.g., 
011 ) allow one to identify the long-distance physics of the loop model with that of the 
sigma model with a topological angle 9 = tt. The loop model sits at a point of hrst- 
order phase transition for n > 2, and is critical for n G [—2, 2]. It is widely believed that the 
sigma model enjoys the same properties, and in fact, the loop model was historically studied 
in part because of the correspondence with Q. 

Another, maybe more physical, way to see the emergence of the U(n) symmetry (see 
[9, 10 for closely related work in both three and two dimensions) is to consider directly the 
Euclidean version, and observe that the partition function can be calculated by introducing 
n-dimensional complex vectors z that live on the edgesjf] with interactions that match the 
geometrical dehnition. This is done by associating, for instance, to the vertex 



a term 


-.-s 


p (4 ■ • ^ d ) + (1 - p )(4 ■ • ^ b ) 


( 3 ) 


where the vectors Tand , living respectively in the fundamental and conjugate fundamental 
representations, obey the normalisation |Tp = = 1, and the bracket represents the 

contribution to the Boltzmann weight from the corresponding vertex (note that there is no 
term without a vector z, since the model is completely packed). The real parameter p can 
be interpreted as the probability of taking the hrst diagram in the expansion of (|^. The 
partition function 

^ oc n J 

edges e ^ 


can be expanded by picking either of the two terms in each vertex contribution, giving rise 
immediately to the loop model with weight n per loop, since the vectors Thave n components; 
the contribution to each vertex is actually proportional to that of (j^, with X = Note 

f The fundamental representation n should not be confused with the loop weight n G C. It is mathematically 
well-defined only for n € N, in which case the dimension of n is n. However, most results throughout this 
paper make sense by analytic continuation for arbitrary n G K. 

§ One may naively wonder why a complex vector z is needed for the loop model. The point is that we need 
the interaction to decompose into two diagrams only, and this is what happens in n(8)h in U(n). In contrast, 
in 0(n), we have three diagrams. This means that an interaction where we would replace the z- z' by n • n' 
could not be interpreted unambiguously in terms of loops. 
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that every Zf. occurs twice, since each edge is shared by a pair of vertices, and the contribution 
along a loop of length k is of the form 




(4 


'^ei 


)• 


This implies that the phase of z^, is not a physical degree of freedom, showing that the target 
is U(n)/U(n — 1) x U(l), that is, the projective space CP"“^. We shall see below how the 
interaction term of the sigma model can be obtained by taking the continuum limit of terms 
such as (|^. To see the origin of the topological term in this picture 11 
in (see footnote 3 of chapter 3). 


one may argue as 


Various modihcations of this model can be imagined, depending in particular on what 
happens to the symmetry. Of special interest is the loop model where crossings are now 
allowed at the vertices and given some Boltzmann weights w, while the loop fugacity 
remains equal to n. Adding such an interaction breaks the U(n) symmetry down to 0(?7,), 
since the crossing ‘split’ does not respect the lattice orientation. The transfer matrix can 
then be written in terms of the projectors on the two independent representations that occur 
in the tensor product of the vector 0(?7,) representation with itself. The universality class is 
modihed, and has been shown to correspond to the low-temperature Goldstone phase of 
the 0(n) model, which is described by the weak-coupling hxed point of the 0(?7,)/0(n — 1) 
sigma model. 

Another familiar modihcation consists in diluting the loop model to allow edges that 
are not covered by loops [^. In general, this dilution is made in such a way that the U(n) 
symmetry is broken, and the loop trajectories do not respect the orientation of the lattice. 
The remaining symmetry is again only 0(n) |^. While the interaction is not the most 
general allowed by this symmetry (since no crossings occur), it is known that crossings does 
not change further the universality class, which is generically 0{n) criticality. 

Meanwhile, a modihcation preserving U(?7,) symmetry can be obtained by allowing next- 
nearest neighbour lines to cross, via the permutation generator Vi^i+ 2 - This is most elegantly 


studied on the triangular lattice, and gives interesting results in particular for n = 0 16 


For other recent developments about loop models, see 10 


Another U(n)-preserving modihcation that has not been studied much consists in 
diluting the loops while preserving the orientation of the underlying lattice (see hgure [^. 
This produces loops which never go straight, and that may or may not meet at vertices, with 
the allowed conhgurations being: 


X,X 





X 

X 


( 5 ) 


The Boltzmann weights are here K per monomer and r per vertex encounter. Our purpose 
in the following will be to study the critical properties of this model, which we shall call the 
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Figure 2: Configuration of the oriented, dilute loop model. 


‘dilute oriented loop model’. In all the rest of this paper, we will put a tilde symbol (~) on 
top of quantities refering to the usual (completely packed) TL model, in order to distinguish 
them from quantities in the dilute oriented model. 

The paper is organised as follows. In sectionj^we discuss the transfer matrix of the dilute 
oriented model in terms of its symmetries, the modules on which it acts, and its centraliser, 
both with free and periodic boundary conditions in the transverse direction. We use to this 
end a mixture of combinatorial and algebraic techniques. In section]^ we map the model onto 
a held theory. We discuss the low-energy limit of this theory and the corresponding critical 
exponents. The phase diagram will be shown to contain a critical line where, remarkably, 
the symmetry is that of a dilute 0(2n) model (although the loop weight is just n). This 
line terminates in a point with SU(n -|- 1) symmetry that is the dilute counterpart of the 
completely packed model. The case n = 1 is special, and we shall show that it contains 
the entire critical line of the six-vertex model, along which the exponents vary continuously. 
While several types of open boundary conditions can be considered, we mostly focus our 
numerical study on the periodic case. However, the conclusions we draw are argued to be 
generic. Finally, section contains the discussion and a few concluding remarks. 


2. Transfer matrix and symmetries 


2.1. The algebra 

Before discussing the possible universality classes of this model, it is interesting to study 
algebraic aspects of its transfer matrix description. This description naturally involves a 


‘dilute’ version of the Temperley-Lieb algebra, but not what is usually called 17,18 the 


‘dilute Temperley-Lieb algebra’, since in our model loop strands are never allowed to go 
straight at vertices. The implication for our dilute oriented loop model is that loop strands 
on odd sites can only be contracted with strands on even sites, just like in the completely 
packed case. 
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The generators of the ‘dilute oriented Temperley-Lieb algebra’ are the following 


p ' — III* 

III I 

123 i 


\J 


I • • • I I 


N 


p — III***! _ I***!! 

■ ■ ■ 1 1 ■ ■ 


( 6 ) 


123 ! Af 

where the dashed lines act as the identity, namely 

; = I + I 

Note that the dilute version of the usual TL generator [cf. ([^] is then obtained as 

+ _ ■ ■ ■ ' \j ' 

P . P ' P III***! ^ l•••ll 

— I I I I r\ i ■ ■ - 


( 7 ) 


( 8 ) 


123 


N 


As for the TL algebra the multiplication relations between the generators follow from 
their geometrical representation. Here are a few sample relations: 

e-Ci =ner, 

6^6+ =ne+, 

e+e+ = 0. (9) 


2,-2. The U(n) symmetry 

The U(n) symmetry (for n G N) of the model in the completely packed was probably hrst 
mentioned by Affleck [^. Associating particles with n possible colours a = 1,..., n on even 
sites, and holes with n colours on odd sites, the Temperley-Lieb generators act on a pair 
of neighbours as {eYai = Sabdcd, and the relation = ne arises simply because of the n 
colours that can propagate along the loops. The choice of particles (resp. holes) on even 
(resp. odd) sites corresponds algebraically to taking the fundamental representation n (resp. 
anti-fundamental h) from the U(?7,) point of view. In the present dilute oriented case, each 
edge now carries the direct sum of the trivial representation, denoted i and of dimension 1, 
and the fundamental (resp. the anti-fundamental). The space on which the transfer matrix 
acts is thus, for integer loop weight n, 

TL = [{i + n) {i + n)]^^ (10) 

The generators can easily be written in coordinates if we associate with the representation 
1 an extra label 0. Introduce now Greek symbols a = 0,l,...,nto describe all the states in 
the trivial and fundamental representations. We have then 

= (7 ~ 5ao)(l “ <^,9o)(l “ <^ 7 o)(l ~ dso)Sai3S^s , 

= <^o0^/3o(l — <^7o)(l — 5<5 o)<575 , 
y VajS = (7 “ '^ao)(l ~ ^/3o)^70^<50<^a/3 • 


( 11 ) 
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Since the e*, ef have non-trivial action only between the two singlets that appear in the tensor 
product (i + n) ® (i + h) (or the product with n, h switched), they necessarily commute with 
the (trivial) action of U(?7,). 


2.3. Combinatorics and modules 


2.3.1. Reminders about the completely packed case. In this section and the few following we 
will consider the case where the ‘space’ direction for the transfer matrix is an open segment, 
corresponding to open boundary conditions for the spin chain. The corresponding geometry 
is a square lattice oriented diagonally with free boundary conditions in the space direction, 
as illustrated in hgure 

For n = q + with q generic (not a root of unity), the representation theory of 


the TL algebra is well known 19 to be semi-simple, with all simple modules corresponding 


to the so-called standard modules Vj. These are indexed by the (even, if the number of 
sites N is even) number j of ‘through-lines’ (or ‘strings’), which are not allowed to intersect 
any arc{j]] In both the transfer matrix and Hamiltonian pictures, the number of through- 
lines is semi-conserved, in the sense that it can only be lowered under the action of the 
TL generators (in other words, the transfer matrix and Hamiltonian have a block-triangular 
structure). Further imposing that pairs of through-lines cannot be contracted with one 
another therefore amounts to considering through-lines as conserved (or, in other words, to 
consider only the block-diagonal part of the transfer matrix or Hamiltonian). 


2.3.2. Standard modules of the dilute oriented case. Through-lines can be dehned in the 
same way in the dilute oriented model, however they are not the only conserved quantities. 
Denoting the parity of edges within any row as eoeo... from left to right, where e stands 
for even, and o stands for odd, let (resp. N^) be the number of even (resp. odd) edges 
covered by the loops within a row. Then 

g=^(iVe-iVo) (12) 

is conserved by the row-to-row transfer matrix. 

The standard modules Vj for this problem are not only indexed by the number j of 
through-lines, but also by their parity. When, say, an even through-line is joined to an arc, 
it must necessarily join to the “odd” end of the arc (on the nearest-neighbour site), and 
hence come out at the other “even” end of the arc. Therefore the parity of each individual 
through-line is conserved by the transfer matrix. The hrst few standard modules are therefore 
labelled as follows 

Vo Vi(e) Vi(o) V 2 (eo) V 2 (ee) V 2 (oe) V 2 (oo) 

We now compute the dimensions of these standard modules by combinatorial means. 

II Note that our convention slightly differs with the usual one in the TL litterature, where the notation Vj 
corresponds to the module with 2j through-lines. 



Dilute oriented loop models 


8 


2.3.3. Ground state sector (no through-lines). Consider first the sector without through- 
lines, namely the standard module Vq. It consists of dilute arc conhgurations, with the 
crucial constraint that each arc connects two sites of opposite parities. This implies that the 
number of empty sites inside any arc A —not counting the sites inside other arcs B contained 
within A —must be even. 

Let t be the weight per site, and let /e(t) denote the generating function for valid 
arc conhgurations on an even number of sites. We similarly dehne fo{t) as the generating 
function for arc conhgurations on a odd number of sites. By examining the possibilities for 
the leftmost site we obtain the functional equations 


/e(t) = 1 +f/o(t) +tVe(^)^ 

fo{t) = tfeit) + t‘^feit)foit) . 


(13) 

(14) 


For instance, the three terms on the right-hand side of (13) correspond to the leftmost site 
being non-existent, empty, or supporting an arc. In the latter case, the arc connects an even 
and an odd site, so the sites inside the arc and those following it are independent (because 
of the self-avoidance) and both described by feit). 

The regular solution for fe(t) reads 




2(-2)V3 ^ (-2)2/3c^V3 


OJ 


1/3 




u 


= - - 4 . 


(15) 

(16) 


The corresponding series expansion is 


fe{t) = 

L=0 


Ql 


1 f3L + 1/ 
L+l\ L )' 


( 17 ) 

(18) 


meaning that ai is the dimension of the standard module Vo for the model dehned on an 
even number of sites 2L. For instance, 02 = 7 and the seven arc conhgurations on four sites 
can be written: 


' \J ' \J ‘ \J ' ' \J \J ViiV 


We note in particular that the two conhgurations • and vj-/ • are forbidden by the 
orientation constraint (the arcs connect sites with the same parity). Those states would 
however be included in the standard module of the usual (not oriented) dilute TL algebra 


17,18 with 0(n) symmetry. 

The similar result for f^it) leads to 


foit) = 


2L-1 


(19) 


b,. = 


L=1 
1 


2L + 1 


3L 

L 


( 20 ) 
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so that now hi is the dimension of Vq on an odd number of sites 2L — 1. For instance, 62 = 12 
and the twelve arc conhgurations on hve sites can be written: 


. \J 'U* U ’ \J ’ ' 

■ ‘ \J \J \J ’ \J \J \J ’ • Viiy' • 


2.3.4- With through-lines. In the above treatment of Vo we have considered systems of 
both even and odd length, because we shall need both generating functions /e(f) and /o(t) 
in the following. For the case with j > 0 through-lines, however, we limit the discussion to 
systems containing an even number of sites 2 L, since this is what is needed to reproduce the 
geometry of hgure 

The generating function for the number of states in Vj{- ■ ■) is the product of the following 


factors: 


• fe for each pair of consecutive parity labels that are different; 

• /o for each pair of consecutive parity labels that are equal; 

• fe if the hrst parity label is even, and /o if it is odd; 

• /o if the last parity label is even, and fe if it is odd; 

• V , where j is the total number of through-lines. 

For instance, V 4 (eoee) is associated with the generating function fff^D. In general, denote 
the sequence of parity labels by (/ii/U 2 ■ ■ ■ p^j), where Pi = 1 (resp. pi = — 1 ) if the Fth label 
is even (resp. odd). We dehne 



( 21 ) 


It follows from the above itemised list that the corresponding generating function is 





( 22 ) 


All cases of an even number of through-lines are therefore described by the generating 


function 



2l,2k+2l 


(23) 


for appropriate values of the integers k and After extensive manipulations this expands 
as 



(24) 


where we note that the coefficient of the initial jg unity. The dimension of the 

standard module Vj (for an even number of sites 2L) therefore reads 



for j = 2k + 21 even, 


(25) 
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where k = ^ = t) for j = 0; and for any even j > 2: 


k= -{j + - 1), 

^ = -(j + £ — /ii + /ij + 1). 


(26) 


In particular, for the case of zero through-lines (A; = £ = 0) we recover from the general 
result (25) the expression (18) for the dimension ol of Vq. 

Similarly, all cases of an odd number of through-lines are described by the generating 
function 

(27) 

with again fc, £ G N. Similar computations lead to the expansion 


^2(k+2i+i) 


4A; -|- 6£ -|- 5 


L=k+2i+2 


L-k-2i-l\L-k-2i-2 


3L + k + l 


2L 


(28) 


where we have singled out the hrst term, where otherwise the unit coefficient can be found 
by analytic continuation of the general expression. The corresponding dimensions therefore 
read 


dk,l — 


4:k+6£+5 ( SL+fc+l \ 
L-k-2e-l \L-k-2i-2) 


a L = k + 2i + l 
ii L>k + 2i+ 2 


for j = 2A; -I- 2£ -I- 1 odd, 


(29) 


where now (26) must be replaced, for any odd j > 1, by 


k = - e + fii - fij - 1), 

£ = -{j + e — p .1 + fij — 1). 


(30) 


2.4- The centraliser 

We give two different constructions of the centraliser and compute its dimension within each 
standard module. The hrst construction is combinatorial and based on the resolution of the 
Markov trace in terms of traces over standard modules. The second construction is algebraic 
and based on oscillator representations of the generators. 


2.4-1. Markov trace construction. In the usual dilute TL case, where the lattice orientation 
is not respected, one can hud the degeneracy Dj of each standard module Vj by inversion of 
the sum rule 

2L 

= (f+i)"y (31) 

j=0 

where dj is the dimension of Vj, and £ the weight of a non-contractible loop. We have (n-fl) 
states on each site, since it can be occupied by any of the n loop colours, or be empty. We 
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shall denote by n (resp. €) the weight of a contractible (resp. non-contractible) loop. It can 
be seen that Dj is a polynomial of degree j in i, 

D, = UjW2 ), (32) 


where Uj{x) is the j’th order Chebyshev polynomial of the second kind. Note in particular 
that Dj does not depend on n. 


The sum rule (31) is insufficient for the oriented dilute TL model investigated in this 


paper, since the number of sectors (standard modules) grows exponentially, not linearly, with 
the size 2L of the system. However, the sum rule is just a special case of the more general 
decomposition of the Markov trace of any element w of the algebra: 


Mtrw = Dj 
vp-) 


)tiv,(...)«>, 


(33) 


where (■ ■ ■) stands for the sequence of j parity labels, and tr denotes the standard (matrix) 
trace. 

We can now hnd all the Dj{- ■ ■) by considering a suitable number of different w. It 
suffices to take words w such that a given subset of the 2L points are linked (by the identity 
operator) to the points immediately above them, while the remaining points are constrained 
to be empty. In the notation of ([^ we have for instance w = 11 I • • • | | | , where the leftmost 
points shown are linked while the rightmost ones are empty. One starts with the word in 
which all points are empty, then the word consisting of one link (residing on a site of even 
or odd parity), then words with two links (of any parities), and so on. The sum rule (31) 
corresponds to the particular case where all points are linked to those above them, i.e., where 
w = I is the identity. 

The matrix traces are evaluated by placing w on top of any basis element b G Vj(- ■ ■) 
and checking if we get back the same basis element, and if so with which weight j^. Such 
a basis element b can be said to be ‘compatible’ with w and hence contributes to trvj(...) 
with the corresponding weight. Since through-lines have a conserved parity, we can give two 
distinct weights, 4 and io, to even and odd non contractible loops. On the left-hand side of 
(33) we get 

(34) 


(4 + 1 )"^ (4 + 1 ) 


Afo 


where (resp. N^) is the number of even (resp. odd) links in w, i.e., the numbers of non- 
contractible loops of each parity in Mtrw. As explained earlier, different through-lines are 
not allowed to be contracted by arcs within the standard module, and the weight of any 
(contractible) loops appearing in the product is n. In fact, with the above choices of w such 
loops actually cannot appear, so all weights are 0 or 1. 

With the above choices of w, it is also easy to see that the basis elements contributing 
to trvj(...)ta are such that all the non-linked points are empty. The linked points must carry 
the number of through-lines with the parity labels specihed by ■ ■ )• Aiiy linked point not 
carrying a through-line can be either empty or carry an arc. Any basis element satisfying 
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these criteria contributes 1 to the trace, so it suffices to count the connectivity states with 
these specihcations. To this end, it is very useful that we already know from the computation 


of the dimensions dj{- ■ ■) = dk/ how to handle the state counting (see section 2.3.4) 


We stress that in these computation all that matters is the appropriate connectivity 
states drawn on the j linked sites. In particular, the total number of points 2L is immaterial, 
and the results for Dj{- ■ ■) are therefore valid for any 2L > j. 

Let us detail one sample computation determining Zl 3 (eeo). We choose any w where 
sites with the chosen labels eeo (e.g., the hrst, third and fourth sites) are linked to those 
above them, and all remaining sites are empty. Thus w = I I | I ■ • •, but it works equally well 
to conduct the computation on just three sites, taking tc = 1 1 1 , provided we take into account 
that the three links have the specihed parities (eeo). We have then trv(,tc = 3 (corresponding 
to the states ' " ', 'si/ or ' \J living on the linked points), trype)"?^ = 3 (corresponding to I ’ 

‘ I ‘ or I VJ on the linked points), trvpo)W = 1 (the only possibility being ‘ ‘ I), trv 2 (eo)'if = 2 
(corresponding to I ' I or ' I l), and similarly trv 2 (ee)'W^ = trv 3 (eeo)'W^ = 1- All other traces are 


zero. We have thus from (33) 


(£e + l)^(f*o T 1) — SDq T 3Z)i(e) + Z)i(o) + 2D2{go) + D2{gg) + D^I^ggo) , (35) 


and since all Dj{- ■ ■) with j < 3 can be found from computations with fewer linked points, 
this eventually determines D^Iggo). 

In this way we hnd 


Do 



= 1 , 


Di\ 

:e) 


= 4 , 


Di\ 

:o) 


— / 


D 2 \ 

(eo) = 

A>2(oe) 

= 44 - 1 , 


D 2 \ 

(ee) 


= D,{g)\ 


D 2 \ 

(00) 


= D,{or, 


Do\ 

(eoe) 


= iX - 24, 


D^\ 

(oeo) 


CD 

0 to 

1 

KD 

d^ 

0 


D3I 

(eoo) = 

= Do{oog) 

= D 2 {go)Di{o 

), 

D^\ 

(eeo) = 

= Do{ogg) 

= D 2 {go)Di{g 

), 

-D3I 

(eee) 


= Ill(e)^ 


D3I 

(000) 


= ^l(o)^ 


D4I 

(eoeo) 

= ^4(0000) 

= iX - 344 

+ 1 , 

D4I 

(eooe) 

= D 2 {Gof . 




From these examples the general result can now be inferred. Whenever the pattern of parity 
labels is alternating (eoeo..., or oeoe...) the result is that of the usual dilute TL model, 
with the obvious replacements i ^ ig or Iq. More precisely, for even j, each monomial has 
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an equal number of and factors, whereas for odd j each monomial has one excess factor 
of £e (resp. £o) for the alternating pattern starting and ending with an ‘e’ (resp. ‘o’). 

Whenever the pattern of parity labels is not alternating, there is a factorisation onto 
contiguous alternating sub-patterns. For example we will have 


DQ^eoeeoo) = D3{eoe)D2{eo)Di{o) 


(37) 


This factorisation is obviously unique (assuming, of course, that each sub-pattern is of 
maximal length). 

We stated in the beginning of this section that the Dj(---) cannot be inferred by 


applying the sumrule (31) alone. However, now that the result has been worked out, it 
is a non-trivial check of the expressions for both dj{-■ ■) and to verify that the 

sumrule is indeed satished. In the case of the dilute oriented model, the sum on the left- 
hand side should carry over both j and the sector labels. For a given number of sites 2L 
we should only sum over sector labels that can be realised from sub-sequences of the basic 
pattern eoeo • 


. Using (25) we hnd, for instance when L = 2, 


7 Z)o(g) T 5 Z)i(e) -|- 5 Z)i(o) -|- 6D2{go) Z)2(oe) -|- D2{gg) 7)2(00) -|- 

Z)3(eoe) -F 7)3(000) -F 7)3(000) -F 7)3(000) -F 7)4(0000) = (7e + 1)^(4 + 1)^ , (38) 


where (36) has been used in the last equality. Note in particular that a term like 7 ) 3 ( 000 ) 
does not appear, because it cannot be formed as a sub-sequence of ‘eoeo’. 


2.4-2. Algebraic construction. We consider the system with free (open) boundary 
conditions, and label the sites i = 0,..., 2L — 1. Recall that we associate an (n -|- 1)- 
dimensional complex vector space U, which is the direct sum 1 © n of the trivial and the 
fundamental for even i (resp. trivial and dual fundamental, i©h, for odd i) representations of 
U(n). It is convenient to use an oscillator representation, and introduce the operators h\^ 
for i even, hia, Vf for i odd, with commutation relations [7“, 7^^] = (with a,b = 1,..., n), 
and similarly for i odd. 

The annihilation operators 7“, hia destroy the singlet state (which can be considered as 
the oscillator vacuum), the daggers indicate the adjoint, and the spaces U are dehned by 
the constraints 

h\ah1 = 0 or 1 [i even), 

b'fhia = 0 or 1 {i odd) (39) 

of having at most one boson per site (we use the summation convention for repeated indices 
of the same type as a). We dehne the generators of U(n) (or in fact of the Lie algebra gU) 
acting in the spaces U by = b\j)\ for i even, = —l\bia for i odd, and the commutation 
relations among the J^’s (for each i) are ©independent. Hence the global gl„ algebra, dehned 
by its generators Jiai acts in the tensor product V = <Sif=Q^Vi. Note that the U(l) 

subalgebra of gl„ generated by does not act trivially on the chain (as it counts the number 
of bosons), in contrast with the completely packed case studied in |^. 
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The Temperley-Lieb generator is the usual ‘Heisenberg coupling’ of magnetism, and can 
be written in terms of the bosonic generators as 


bf+Aj)%+i,h, i even, 
Abl+i,abt+Ab, i odd. 


(40) 


The Cj’s are Hermitian, e| = e*. Acting in the constrained space V, they satisfy the 
relations 

= net, 

dej =ejei {j^i,i±l). (41) 

While Ci acts as e* = 0 whenever the number of bosons on one of the two neighbouring spaces 
is zero, we can also introduce the generators 


and 


b%+i^b, i even, 
i odd. 




i even, 
i odd. 


(42) 


(43) 


which act as required in the constrained spaces. The are Hermitian conjugate, and 
U(n)-invariant. 

To construct the commutant algebra explicitly, we introduce the operators (for k < 2L) 


/T-aia2...afc 

bib2...b^. 



ja\ ja2 _ _ _ TOfe 
^ilbi^i2b2 


0<ii<*2<"-<*fc<2L—1 


(44) 


(for /c = 0, we dehne = 1, and for /c = 1, as dehned earlier). We moreover impose 

linear conditions, that the contraction of one of the indices a with a neighbouring index b 
[i.e. of ai with 6;+i (resp., &z_i), for Z = 1, 2, ..., fc — 1 (resp., Z = 2, ..., Zc)] is zero. This 
gives us a basis set that is ‘traceless’ in this sense. For example, for k = 2, we have 


jCLlCL2 _ rraia2 

^bib2 ~ '-^h. 


'bib2 


_ 'Ta,a 2 j:ai 

"62 


n 



(45) 


and there are {n^ — 1)^ independent such operators. In general, there are {D^Y of them, 
where 

= [3 + ( 46 ) 

and where [n]q = + ... + = (g^ — q~'^)/{q — q~Y is the g-deformation of 

any integer n. 

The exact forms are 

( p» p rj-\axa2...ak 
W )bib2...bk ’ 


ja\a2...ak 

'^bib2...bk 


( 47 ) 
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where P* (resp. P,) is the (Jones-Wenzl) projection operator to the ‘traceless’ sector on the 
vector space indexed by (oi, 62 , • • •) [resp., ( 61 , 02 ,...,which can be constructed recursively 
using the TLfc(q') algebra in these spaces (see for a review). 

In the ordinary (completely packed) case, it is known that the generate the 

centraliser of the TL algebra in the sector with k through-lines. The total dimension of the 
centraliser An{2L) is then 

L 

dim An{‘2L) = . (48) 

i=o 


Note that dimensions D 2 j are independent of L, so the total dimension of the centraliser is 
determined by the upper bound in the sum, j^ax = L. The centraliser is simple for n an 
integer n >2, and the D 2 j are the dimensions of its irreducible representations. 

The centraliser in the dilute oriented case is then obtained by acting with the generators 
through lines that can be contracted by the algebra. This set of course 
depends not only on the number of through lines, but now also on the associated pattern, 
since only lines on odd and even sites can be contracted. For instance, for the sector ‘oeoe’ the 
associated representation of the centraliser has dimension iJ 4 = ZJ 4 (oeoe), while for ‘oeeo’ it 
has dimension D 2 and for ‘oeee’ it has dimension D 2 D 1 . This factorisation phenomenon has 
already been discussed above; see (37) in particular. In general the dimension is a product 
of Dki with '^ki = j, the total number of through lines. We will denote the centraliser in 
this case by An{2L). Its dimension is given by an expression analogous to (48), except that 
the sum has to carry over both j and the associated sector labels that arise as sub-sequences 
of the basic pattern eoeo ■ ■ ■. For instance we hnd, in analogy with (38), that for L = 2 


dim^„(4) = Dl + 2Dj + 2Dl + 2{DiDiY + 2Dl -F 2{DiD2) 
= n® — 2 n® J- 3 n"^ J- 2 n^ -\- 4 . 


+ Di 


(49) 


2.5. The limit {K,t) —)■ (|,2) and the SU(n-|- 1) point 

Something remarkable must happen at the particular point where the parameters K and 
r introduced in equation ^ are {K,t) = (|, 2 ), where the model exhibits an extended 
SU(n-|-l) symmetry. This is most easily seen starting from the completely packed model with 
loop fugacity n-l-1, which obviously has an SU(n-l-l) symmetry. Imagine now colouring each 
loop in two possible ways: either black (corresponding to n loop components) or ‘transparent’ 
(corresponding to the last component) (^. The resulting configurations have a fugacity n 
per black loop, while the fugacity of transparent loops is trivial, so that these need not 
be counted in the partition function and can hence be considered invisible indeed. As a 
consequence of this transformation the black loops appear diluted. This diluted model 
coincides exactly with the model we are discussing in this paper when iF = 4 and r = 2 
(consider twice the right-hand side of (|^, bearing in mind that the hrst ‘empty’ diagram 
corresponds now to the two different splits of the transparent loops). 
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In algebraic terms, what must happen is that the centraliser of the algebra, ^„(2L), 
must become bigger at this special value of the coupling r, and coincide then with An+i{‘2,L). 
In other words, representations of An{‘2^L) combine to form representations of An+i{‘^L). This 
must coincide, of course, with the appearance of extra degeneracies in the transfer matrix 
spectrum, such that the decomposition of the Hilbert space of total dimension (n + 1)^'^ 
can be accomplished using both the dilute oriented TL algebra with parameter n and the 
ordinary TL algebra with parameter n+1, as well as their commutants v4.n(2L), and An+i{‘2,L) 
respectively. 

For the ordinary TL algebra, we decompose [(i © n) ® (i © h)]®'^ as 

L 

{n + lf^ = ^b 2 jd 2 j, (50) 

j=0 

where the Dj are the dimensions of the centraliser 

Dj = [j + l]g) [2]g = n + l = n. (51) 

We have for instance Di = n, D 2 = = rri — 2n and D 4 = — 3rP + 1. For the 

dilute oriented TL algebra and its centraliser, we will have a different and more complicated 
formula, which we will analyze explicitly for L = 2. We hrst introduce the notation 


= li + 1] . [2] = n , 


(52) 


and rewrite (38) as 

{n + 1 )^ = do T [2](di(e) + di(o)) + [ 3 ](d 2 (Go) + d 2 (oe)) + [2]^((i2(ee) + ^2(00)) 
+ [4](d3(oeo) + d 3 (eoe)) + [2][3](d3(ooe) + d 3 (oee)) + [5]d4(oeoe) .(53) 

This reads explicitly 

(n + 1)^ = 7 + (5 + 5)u + (6 + l){n^ - 1) + (1 + l)n^ 

+ (1 + l)(n^ — 2n) + (1 + 1)(?T'^ — n) + (n^ — 3n^ + 1). (54) 

For the ordinary TL we would have instead 

{n + lf = n^ = 2 + 3{n^ - 1) + + 1) . (55) 


For this to be possible, it is necessary that eigenvalues from the dilute oriented model acquire 
extra degeneracies at the point {K,t) = (|,2) so that the dimensions of the commutant 
become larger. Still in the case L = 2 we hnd that 

3(^2 _i) =3 + 6n + 3(n2-l) 

n'^ — 3n^ + 1 = 2 + 4U + 4(n^ — 1) + + 2(n^ — 2n) + 

2(n^ — n) + ijb — 3n^ + 1), 


( 56 ) 
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so we see, for instance, that in the representation of dimension do = 7 ol the constrained 
TL algebra, 2 eigenvalnes will remain non degenerate, 3 will become degenerate with some 
of those in the sector with one and two non contractible legs, and 2 will become degenerate 
with some of those in all the other sectors. This is indeed what we check nnmerically for 
generic valnes of n, and can be snmmed np as follows: 


do —7 

di{e) + di{o) = 10 

d2(eo) + d 2 {oe) = 7 

d 2 (ee) + ^2(00) = 2 

doi^eoe) + ^ 3 ( 000 ) = 2 

do^ooe) + do^eeo) = 2 

di{eoeo) = 1 


(do — 2) X + 

i I 


{d 2 — 3) X 

1 

2 
1 


+ (^4 — 1) X 

2 
4 
4 
2 
2 
2 
1 


where the numbers appearing along the vertical bars correspond to the multiplicities with 
which each dense module enters the decomposition of the dilute modules. Finding the 
corresponding branching rules for general L is an interesting algebraic exercise, which is 
however beyond the scope of this paper. 


2.6. The periodic case 

A similar procedure could be repeated in the case of a lattice with periodic boundary 
conditions in the horizontal direction. 

As in the dense case, some care has to be taken concerning the treatment of loops 
wrapping around the space-like direction. Namely, in the sector with zero through-lines, an 
arc between two given points may or may not intersect the periodic boundary condition (a 
vertical ‘seam’). These two possibilities (which we can depict as and ^ v in terms of the 
usual diagrams) can be identihed—formally, by taking an algebra quotient—or considered 
different, leading in the latter case to additional multiplicities. In the following we will focus 
on the former (quotient) case. In the sectors with j through-lines, one must also attribute a 
dehnite momentum to their winding around the seam. Specihcally, each through-line picks 
up a phase ip = (with m = 0,1,..., j — 1) when it crosses the seam towards the right 
(and —(p when it crosses the seam towards the left). Finally, one must notice that in this 
periodic case modules whose parity labels differ by a cyclic permutations—for instance, the 
modules ‘eo’ and ‘oe’—are to be identihed. 

We will not go into any computational detail here, but simply observe that at the 
SU(n -|- 1) point {K,t) = (^,2) the spectrum in each module can similarly be decomposed 
in terms of that of the periodic dense model. For L = 2, for instance, we hnd the following: 
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do =7 

di(e)(°)+di(o)(°) = 12 

d2(oe)(°^ = 8 

d2(oe)(^) = 8 

^2(00)^°^ + c? 2 (ee)^°) = 2 

d2(oo)W+d2(ee)(") = 2 

^ 3 ( 000 )^°^ + d 3 (eoe)^°^ = 4 

^ 3 ( 060 )^^^ + d 3 (eoe)(^^ = 4 

d4(eoeo)(°) = 1 

d4(eoeo)(l^ = 1 

d4(eoeo)(’^^ = 1 


(do = 2)x 
1 


+ (rff=4)x 
1 
2 
1 


+ (4"^=4)x 


1 


+ (^4 ^ — 1) X 

1 

2 
2 
2 
2 


2 

2 

1 


1 


+(4^|=i)x 


2 

2 

2 

2 

2 

2 

1 


where for each module with through-lines we have indicated between parentheses the value 
of the corresponding phase 99 , as dehned above. 


3. Low-energy limit and critical points 

3.1. Mapping onto afield theory 

In order to understand what kind of critical point one may observe in this model, it is useful 
to start by considering in more detail the map of the completely packed theory onto the 
CP"“^ model. We thus go back to (|^, introduce a continuous variable z{x,y), and expand 
this variable around the vertex. For instance we write 

a? 

za = z + a(-d^z + dyZ) + y dyZ- 2d^dyZ} + ... , (57) 

where a is the cutoff and we restricted to second order. A painful but straightforward 
calculation gives 

{^A • ^b){4; ■ zd) = 1- 4a^ [{z- dyZ^f + {dyZ.dyZ^)] + ... (58) 

and similarly for the other term up to the exchange x ^ y. At the isotropic point p = 1/2 
to which we will restrict, we see that the Boltzmann weight e""^ can be re-exponentiated, 
giving 

S'oca^ ^ {z-dy^fi + {dyZ-dy^)^ j dxdy [{z ■ + {d^z ■ dy^)] , (59) 

vertices 

where fi = x,y and the sum over both directions is implicit. This action can be checked to 
coincide with the standard action for the model ms 

S ^ J dxdy{D^zt^{D^zt. (60) 

^ The bare coupling is of order 0(1) and entirely determined by the model, which has no free parameters 
(at the isotropic point p = 4). Of course, this coupling flows under renormalisation. 














Dilute oriented loop models 
where the covariant derivative is 


19 


= ^ ■ d^z. (61) 

This action is invariant under the gauge transformations z{x,y) —)■ ?/). 

We now consider what happens for the dilute oriented model. We must add to the 
Boltzmann weight expansion (|^ a new kind of term, which corresponds to the possibility of 
having no loop or only one arch of a loop. We thus take now 


e-^ = 


+/3 


1 +a • zb){^c • ^d) + ■ zd){^c ' ^b) 

(4 • ^b) + (4 ■ ^d) + (4 ■ ^d) + (4;-Zb) } • • • 

The same calculation as before gives now 

(4 ■ ^b) + (4 • = 2 - ia^d^z- , 

so we see that, if we re-exponentiate in the action, we will get generally 

S «. (ixdy [a(D^z)t(£)pZ) +. 


(62) 


(63) 


(64) 


This is a model with two parameters, in agreement with the fact that we have two parameters 
t,K on the lattice. The important fact is that the new term breaks gauge 

invariance. 

The situation we are encountering has in fact a long history in the area of sigma models, 
in particular in what is called ‘isovector-isotensor’ models. For instance, in 24 a model with 
the following lattice hamiltonian is studied 


fdy^x 


I3t 


' 2 


(J',, 


(65) 


" = -E 

i^y) 

where a is an n-component real spin obeying = 1. The case = 0 is the usual ‘0(n) 

model’ (or, more precisely, the model with target 5'”'“^ = 0 (n)/ 0 (n — 1 )), while the case 
/Sy = 0 is the model. This latter case has an extra local (gauge) symmetry, since 

we can change the sign of a at each point independently]^ 

Unfortunately there are not many results available on the phase diagram of similar 
deformations for the model. In the RP”“^ case, it seems that the physics of the 

model with both /9y,/9r 7 ^ 0 is the same as that of the ordinary 0(n) model 24 . Since 

RP"-^ 


RP 


n-l 


= S'"' ^IZ 2 ., we see that mixing of the isovector term renders the target bigger, 
—)■ By analogy, we expect that in the case of complex vectors we have 


CP"-^ = SU(n)/SU(n - 1) 0 U(l) ^ SU(n)/SU(n - 1) = R 


2n-l 


( 66 ) 


that is, one should expect to observe—at least in some regions of the phase diagram—the 
physics of the 0 ( 2 n) model. 

+ Note that in the usual 0(n) model on the hexagonal lattice, once one has restricted to closed planar loops, 
every spin S occurs twice, and one can also change its sign at will. But this 0(n) model is not the same as 
the RP”~^ model because the action is not gauge invariant. The same occurs in our case: if we keep only 
the loop diagrams, the phases of the z still cancel out. But the continuous action is not gauge invariant. 
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1/2 


Figure 3: Qualitative phase diagram of the oriented dilute loop model for some hxed value 
of u < 1. For small monomer fugacity K the system is in a massive phase with zero average 
density of monomer in the thermodynamic limit. For large values of K one enters the so- 
called ‘dense’ phase, which is critical and described by the dense 0(n) universality class. 
The two phases are separated for low r by a critical line (thick red curve), and for large r by 
a hrst-order line (double dashed curve), in between which sits the SU(n -|- 1) point (purple 
dot). 


3.2. Continuum limit at the SU(n -|- 1) point 

The phase diagram of the dilute oriented loop model in the {K, r) plane for any — 1 < n < 1 
is in fact quite simple. Its shape is represented schematically in hgure and it exhibits in 
particular a single critical line separating a dense and a dilute phasej^ 

Before examining this critical line in detail, we focus in this section on the SU(n -|- 1) 
point {K,t) = (|,2); see section 2.5 It is indeed straightforward to understand the critical 
content of the dilute oriented model at that point in terms of the usual SU(n -|- 1) dense 
critical model. Let us hrst recall a few facts about the latter. Parameterising 


n = n + 1 = —2 cos ttq , with 0 < < 1, 


(67) 


the central charge—which is related to the hnite-size scaling of the ground state (i.e., the 
lowest-energy eigenstate of the Hamiltonian in the sector with zero through-lines)—reads 


c = 1 — 6 


(1 - 
9 


( 68 ) 


and the critical exponents—which are related to the hnite-size scaling of the gaps between 
excited eigenlevels and the ground state—can be expressed in terms of a Kac-like formula 


hr.s 


{sg-rr-{l-gf 

^9 


(69) 


* In the following we consider that the continuous variable n takes some fixed, given value, and discuss the 
phase diagram in the {K, t) plane that intersects this value of n in the full (n, iti, r) space. In particular, we 
shall refer to the SU(n + 1) “point” (although it is rather a “line”, as a function of n). 
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In particular the leading exponent in the sector with i = j through-lines (the so-called 
watermelon exponent Xi) reads, in the open case [^[2^ 


xe — — 


4:9 


(70) 


and in the periodic case 27 


Xi = 2hg e = 


~9 {l-~9f 


2~g 


(71) 


We now go back to the dilute oriented model. Quite generally, we can similarly dehne 
scaling exponents Xe, Xo, Xeo, ■ ■ ■ associated with the scaling of the leading eigenvalue (lowest 
energy) in each sector of j through-lines with prescribed parities||] At the special point 
(77, r = |,2) these exponents should be related to the Xi of the dense model. We have 
checked numerically that the correspondence reads 


Xq — Xq — X 2 1 


X. 


eo 


^oe ^2 
Xqq ^00 ^4 


(72) 


and so on. The general rule is that, given a sequence of parities in the dilute model, 
the corresponding number of through-lines in the completely packed model is obtained by 
hlling in the minimal number of through-lines such that the sequence becomes alternating 
(eoeo- 


eo). This is also in agreement with the decompositions in sections 2.5 and 2.6 


In conclusion, the exponents of the dilute oriented model at the special SU(?7. -|- 1) point 
are the watermelon exponents X 2 j of the associated completely packed model, where the 
(even) number j of through-lines is dictated by parity constraints. 


3.3. Numerical evidence for the 0(2n) universality class 

We now turn back to the study of the critical line below the SU(rr-|-1) point, that is, the red 
line in £gure|^for generic values r < 2. We also shall restrict from now on to the bulk case, 
corresponding to a periodic Hamiltonian. In the boundary case different kinds of critical 
behaviour can be observed, depending on the weight given to boundary loops 28,29 . In 


conjuction with the already quite complex bulk behaviour, which we shall describe below, 
this makes a clear convergence of the critical exponents to any of the natural candidates 
very difficult to assess, at least for the range of system sizes considered here (namely, up to 
2L = 18). Accordingly, we shall leave the examination of the boundary criticality to future 
work. 


Clear evidence for the 0(2n) critical behaviour suggested by the analysis of section 3.1 
is obtained by looking at the central charge, whose numerical estimation as a function of n 


tt For simplicity, we suppress j in this notation, since its value can be inferred by counting the number of 
parity labels. 
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• L=4 

• L=6 

• L=8 

• L=10 

• U12 


Figure 4: Central charge on the dilute critical line (we took here r = 1), measured from 
hnite size scaling of the transfer matrix eigenvalues at sizes L,L + 2 for L = 4, 6, 8,10,12, 
plotted as a function of n. The data shows remarkable agreement with the 0(2n) central 


charge (73), plotted in comparison. 


shows remarkable agreement with that of the 0(2?7,) dilute branch, namely 


c = 1 — 6 




2 n = —2 cos{Tig ), 1 < < 2 . 


(73) 


(see £gure|^. Here and below we place a circumflex (") on top of qnantities referring to the 
dilnte 0(2n) model. We note that values of the central charge for the particular case r = 1 


were already displayed in Table V of 30 ; the relationship with dilnte 0(2n) criticality was 
however not noticed in this reference. 

To proceed further, we recall that the bulk watermelon exponents at the 0(2n) critical 
point read 


xe = 


^9 {i-gf 


2^ 


(74) 


In the following sections, we will give nnmerical evidence for the fact that the exponents 
of the dilute oriented model on its critical line—that is, on the red curve K = Kc{t) in 
hgure 1^ — are indeed given by ([74l) for any valne of r < 2. Here i simply equals the number 


of through-lines, regardless of their parities, and the parity constraints fonnd in the above 
discussion of the SU(n -|- 1) simply disappear on the dilnte critical line. This means that 


(72) is replaced by 


Xe= Xo = Xi 

rf> - - rf> - 'lf*n 

0/00 — — o /2 


(75) 


In spite of the nnmerical problems evoked above, we would expect—by analogy with 


the bulk case—that the exponents with open bonndary conditions are still those (70) of the 











Dilute oriented loop models 


23 



SU(n + 1) point 


dilute critical line 


Figure 5: Evolution of the dilute oriented model’s critical exponents from the SU(n + 1) 
point to the dilute critical line, in the case of periodic boundary conditions. The exponents 
xi and xi are given by (|7l|) and (74) respectively. An analogous diagram can be drawn i 


m 


the open case, with the exponents xi and Xi being now given by (70) and (76) respectively. 


ordinary surface transition, which we rewrite here as 

[(£ + 1 )^- If -( 1-^)2 


Xi = 


4 ^ 


(76) 


with 2 = 1 and no parity constraints. The situation is summarised in hgure 


Let us proceed to check the prediction (75) for the bulk exponents numerically for 


arbitrary values of n < 1 and r < 2. For reasons that will be made clear in the following, it 
is convenient to hrst look at values of n and r far enough from n = 1 and r = 2, respectively. 


3.3.1. First look: negative n, small r. We therefore take, for instance, n = — | and r = 0, 
and compute the critical exponents by exact diagonalisation of the loop transfer matrix, by 
varying K across the critical value K = Kc{t), for system sizes ranging up to 2L = 16. 

As shown in hgure the convergence of the exponent Xe = Xo at the critical value 
K = Kc{t) is in perfect agreement with the 0(2n) value Xi given by ( [74l ). While for the 
exponents x^e = Xoo and Xeo = Xoe the convergence is not as good (see hgures[^and|^, htting 
the hnite-size results to a quadratic function of L~^ leaves little doubt that these exponents 
are once again described by (75). 


3.3.2. Increasing n and r. After these promising conclusions for negative n, we now look 
at larger values of both n and r, for instance n = ^ and r = 1. As shown in hgure where 
the exponent Xg = Xo is studied, the convergence towards the expected Xi is, at best, very 
slow. The analysis is based on the crossings near K k, 0.633; the hgure exhibits another set 
of less neat crossings around K k. 0.69, but their physical relevance can easily be discarded 
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Figure 6: Watermelon exponent = Xq across the n = — | dilute to dense transition, at 
r = 0, obtained from exact diagonalisation of the periodic transfer matrix for system sizes 
2L = 4 (bottom curve on the left of the graph) to 2L = 16 (top curve on the left on the 
graph). We show for comparison the exponent X 2 for the dense 0(n + 1) model (in black), 
and the exponent hi = for the dilute 0(2n) model (in red). Convergence towards the 
latter is almost perfect. 



Figure 7: Watermelon exponent Xoe across the n = — | dilute to dense transition, at r = 0 
(periodic boundary conditions). We show for comparison the exponent X 2 for the dense 
0(n+ 1) model (in black), and the exponent X 2 for the dilute 0(2n) model (in red). Looking 
more closely at the intersections in the left hgure one sees that all curves corresponding to 
different sizes do not intersect at the same location. So in the right panel we study the 
behaviour of the intersection of the exponents for sizes 2L and 2L + 2 as a function of L. 
The blue curves hts the last three sizes {2L = 10,12,14) to a parabolic function, and yields 
an estimate at L —)■ cx) which coincides precisely with X 2 = ^ vo the dilute 0(2?7,) model. 
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n=-0.5, r=0, p=0 



Figure 8: The same as figure but for the exponent Xee- The black (resp. red) lines shows 
Xi for the dense 0(?7, + 1) model (resp. X 2 for the dilute 0{2n) model). The right panel 
provides the intersection of the exponents for sizes 2L and 2L + 2, and the blue curves fits 
those values for 2L = 4, 6,8,10,12 to a parabolic function. The resulting estimate at L —)■ cx) 
coincides precisely with a ;2 = ^ in the dilute 0(2n) model. 


n=0.5, r=1, p=0 



n=0.5, r=1, p=0 



1 

2L 


Figure 9: Watermelon exponent x^ = Xo across the n = ^ dilute to dense transition, at 
r = 1 (periodic boundary conditions). We plot for comparison the exponent X 2 for the dense 
0(?7, + 1) model (black line), and the exponent Xi = | for the dilute 0{2n) model (red line). 
The finite-size analysis in the right panel proceeds in analogy with the preceeding figures. It 
shows that as L —)■ oo the exponent decreases, very slowly, presumably towards the dilute 
0(2n) value. 
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n=1 



Figure 10: Location of the critical point (iFc, t'c) in the n = 1 plane, estimated numerically 
for several values of r (blue dots). We plotted in comparison the line tK"^ = which indeed 
seems to parameterise the line of critical points. 


by examining the central charge. The low quality of the convergence originates, we believe, 
in the fact that for n = 1, the 0(2n = 2) model has, in fact, a line of critical points, along 
which the exponents are expected to vary continuously. Let us now examine this scenario. 


3.4- The special case n = 1 

Our numerical determination of the critical line Kc{t) in the plane n = 1 is shown in hgure 
[To| We already know analytically that it contains the SU(1 + 1) point at {K,t) = (|,2). 
The figure gives convincing evidence that (iFc)^ = ^ for any r G (0,2], and we shall show 
analytically below that this is indeed the true answer. Moreover, we shall show that the 
critical exponents vary continuously with r, and give the exact values for a subset of the 
exponents along the critical line. 

First, we notice that the model for n = 1 can be mapped onto an eight-vertex modejtt] 


by assigning arrows to the edges as explained in figure 11; for convenience we have rotated 
the vertices by 45° in the hgure. 

The last two pictures correspond to the same eight-vertex conhguration, with total 
weight 2K‘^t. There is thus a total of six vertices out of the eight possible that have a 
non-zero weight. 

Any vertex model on a bipartite lattice is invariant under a gauge transformation 
consisting in reversing the arrows on the South and East edges (resp. the North and West 
edges) on the even (resp. odd) sublattice. We now apply this transformation to the eight- 


vertex model in hgure 11 The result is that we recover a staggered six-vertex model, where 


ft This mapping and the subsequent mapping to a six-vertex model were already discussed in 30 , but only 
for the particular case {K,t) = {^D)- 
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1 K K K K K^t K^t 


Figure 11: Mapping from the dilute model in the n = 1 plane onto the eight-vertex model. 
Occupied (resp. empty) horizontal edges correspond to an arrow towards the right (resp. 
left), while occupied (resp. empty) vertical edges correspond to an up-pointing (resp. down¬ 
pointing) arrow. The corresponding Boltzmann weights are also indicated. The two omitted 
eight-vertex conhgurations have zero weight. 


the weights on the even sublattice read (in the standard notation [^) 

{UJI,UJ2, UJ3, uiA, (^5, = (K, K, K, K, 1, 2iFV), 

while those on the odd sublattice are 


(77) 


(78) 


We recover a homogeneous, and hence solvable, six-vertex model if and only if ooi = cu' for 
all i, that is, 2K‘^t = 1. The corresponding anisotropy parameter is then (still in the usual 
notation 1^) 


A = 


+ b^ 


2ab 


= 1 — r . 


(79) 


When r varies between 0 and 2, A varies between 1 and —1, that is, we cover the whole 
critical line of the six-vertex model. For r > 2 the model remains solvable, but since A < — 1 
it belongs to a non-critical phase whose correlation length decreases monotonically with r . 

We note that the gauge transformation preserves the periodic boundary condition, so 
the six-vertex model enjoys purely periodic (not twisted) boundary conditions. The central 
charge in the critical regime r G (0, 2] is therefore that of a free bosonic held, namely c = 1 
all along the line. This is indeed what we hnd numerically. 

The conformal dimensions can be written in the Coulomb gas setup as 


XE.M 


(80) 


where E and M are the so-called electric and magnetic charges (the latter is related to the 
six-vertex magnetisation, M = Sz)- We have here parameterised 


A = cos Tig , 


( 81 ) 
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using now a bar (“) to distinguish the present theory from the two other Coulomb gases used 
so far. 

Let us examine the leading eigenvalues in each sector of the dilute oriented model in 
the different sectors {Sz = 0, ±1,..., ±L) of the six-vertex model. For small sizes, we check 
the following: 

0 —^M = 0,E = 0; 

e = o —)■ M = ±1, = 0 ; 

ee = oo —?• M = ±2, = 0 ; 

eee = ooo —> M = ±3, = 0 ; 

eo = oe —> absent from the six-vertex spectrum; 

eoe = eeo = ooe = ... —)■ absent from the six-vertex spectrum. (82) 

which we interpret as follows. First, it is easy to see from the explicit mapping how each 
non-contractible parity label ‘e’ (resp ‘o’) results in a —1 (resp -|-1) contribution to the total 
value of Sz- The sectors involving parity labels of only one type are in this sense ’highest 
weights’, and are correspondingly found to make up the leading spectrum of the six-vertex 
model. So the states Xq, Xqo, a^ooo, ■ ■ ■ correspond to M = 1, 2, 3,... and E = 0—and similarly 
for Xe, Xee, Xeee, • • •, upon changing the sign of M. Therefore 

_ f . . 

^oo...{l times) 

a result for which we hnd very good numerical evidence (see hgure [I^ , and which recovers in 
particular the known exponents at the SU(1 -|-1) point (which corresponds to (FT, r) = (|, 2), 
so that A = — 1 and namely ^ = 1). 

Conversely, the sectors mixing the two parity labels are found to be absent from the 
six-vertex spectrum, and our interpretation is that the associated operators are non-local 
in the six-vertex formulation. In practice, we observe that the exponents associated with 
purely alternating parity labels are constant all along the critical line, and given by their 
value at the SU(2) point. For instance, we hnd Xeo = Xoe = x^oeo = Xoeoe = 2, as displayed 
in hgure and we conjecture that in general 

Xeoeo...(£ pairs) ^oeoe...(£ pairs) 

As for the remaining exponents, namely those with parity labels mixing even and odd labels 
in a non purely alternating way, we could not produce any analytical formula. From the 
results in hgure the exponents clearly vary along the critical line, however we could not 
formulate a convincing conjecture. 

3.5. The special case n = 0 

The case n = 0 is very special, since formally 0(2n) = 0(n). This means that, starting with 
dilute oriented loop model where loops get vanishing fugacity, we are able to reach the same 
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n=1 



r 


Figure 12: Watermelon exponents of the dilute model along the critical line for n = 1 
(periodic boundary conditions). We have plotted in comparison the expressions (83) for 
^ = 1 and i = 2, with which the exponents with pure parity labels show perfect agreement. 
Meanwhile, the exponents with purely alternating parity labels are constant, namely equal 
to their value at the SU(2) point, all along the critical line. 


universality class as the ordinary dilute model—that is, the model which does not respect 
the lattice orientation. In other words, the universality class on the red line of hgure (|^ is 
the usual dilute 0(0) universality class. Meanwhile, the critical point (violet dot) is in the 


universality class of the theta point as identihed in 31 


4. Conclusion 


The main result of this paper is that, while the ordinary dilute loop model with fugacity n 
per loop has a critical point in the 0(?7,) universality class, the dilute oriented loop model has 
a critical point in the 0(2?7,) universality class. There are, associated with this observation, 
many interesting algebraic as well as phenomenological aspects which deserve further study. 
In particular, note that the dilute oriented model is an special case of a more general loop 
model—namely the usual, non-oriented loop model on the square lattice 15 —that possesses 


and intriguing phase diagram whose features are not all understood at the moment, despite 


some recent progress 32,33 


Of particular importance is the behaviour for n = 0: only in this case are the 0(n) and 
0(2n) universality classes identical, and the underlying lattice orientation mostly irrelevant— 
at least in the continuum limit. In [^, a loop model similar (albeit with more complicated 
rules, involving in particular two loop colours) to the one discussed here was extended from 
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Figure 13: Watermelon exponents a^eoo, x^oee-, Xeoeoo of the dilute oriented model along 
the critical line for n = 1 (periodic boundary conditions), estimated from transfer matrix 
diagonalisation for sizes L = 6, 8,10,12,14,16 (coloured dots). The black dots represent an 
extrapolation to L —>■ oo using a quadratic £t in 
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a definition compatible with the lattice orientation to one that is not, and it was argued that 
the universality class is not modihed upon breaking the lattice orientation. This is not at 
all an obvious result, and most likely holds only for special values of the loop fugacity, just 
like in the model we studied in this paper. This is an aspect we will discuss more elsewhere. 

Meanwhile, the reformulation of the SU(n + 1) model as a dilute oriented model with 
S'U(n) symmetry turns out to be an important step in the study of loop reformulations of 
critical points between universality classes of topological insulators, an aspect which we will 
also explore in more detail elsewhere. 
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